Multiphoton transitions between energy levels in a phase-biased Cooper-pair box 
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We investigated both theoretically and experimentally dynamic features of a phase-biased charge 
qubit consisting of a single-Cooper-pair transistor closed by a superconducting loop. The effective 
inductance of the qubit was probed by a high-quality tank circuit. In the presence of a microwave 
power, with a frequency of the order of the qubit energy level separation, an alteration of the qubit 
inductance was observed. We demonstrate that this effect is caused by the redistribution of the 
qubit level population. The excitation of the qubit by one-, two-, and three-photon processes was 
detected. Quantitative agreement between theory and experimental data was found. 

PACS numbers: 74.50.+r, 85.25.Am, 85.25.Cp 



I. INTRODUCTION 

During the last decade a number of proposals for con- 
structing an artificial quantum two-level system by mak- 
ing use of mesoscopic Josephson junctions were imple- 
mentediiSiiiii. Since it was recognized that these cir- 
cuits might serve as quantum bits (qubits) for quantum 
information devices the field has attracted increased at- 
tention. 

Basically, two kinds of such devices have been devel- 
oped, based on the charge or flux degree of freedom. Here 
we will consider one realization only, which is based on 
the charge degree of freedom. In this quantum system 
two charge states differing by 2e (e is the electron charge) 
are mixed by Josephson tunneling. One example of such 
a device is the single-Cooper-pair transistor - two meso- 
scopic tunnel junctions separated by a small supercon- 
ducting island on which the charge can be induced by an 
external gate voltage^. The relative energy of the states 
is controlled by the gate voltage. 

Since the measurement of a quantum system is a very 
delicate procedure, the readout sensor is a crucial com- 
ponent of any potential quantum computing circuit. In 
order to minimize the exchange of energy between detec- 
tor and qubit the control of the reactive component of 
the output signals has been proposed and implemented^. 
Such kinds of measurements requires the proper design 
of the qubit. For instance a charge qubit can be a con- 
ventional single-Cooper-pair transistor closed by a super- 
conducting loop^. For a certain range of the relationship 
between effective Josephson coupling and charge ener- 
gies Sj/Ec of the transistor's junctions, this device is 
effectively a two-level quantum system with externally 
controlled parameters^'^'^'^^. Moreover, similar to both 
the traditional nonhysteretic RF SQUIDiiii^ and the DC 
SQUID^^, the phase-biased transistor coupled to a high- 
quality radio- frequency tank circuit turns out to be an 
ideal parametric converter of charge and flux signals with 
standard quantum limit of the energy resolution. 



Recently measurements of the energy level separation 
of a superconducting charge qubit were reported. The 
qubit was coupled to a high quality tank circuit^^ or non- 
resonantly to a single mode of the electromagnetic field 
of a superconducting on-chip resonator^^'^^. Multipho- 
ton transitions between energy levels in superconduct- 
ing devices were studied in several article3i2ii2i2Si2ii2S; in 
this work we present both the experimental observation 
and the theoretical description of the multiphoton tran- 
sitions between the ground and the first excited state in 
the phase-biased charge qubit^^ ^^^i-*^^ making use of the 
impedance measurement technique^SiSi. 

We begin in Sec. II with a theoretical description of the 
phase-biased charge qubit (PBCQ) subjected to a time- 
dependent gate voltage or magnetic flux. We calculate 
the population of the upper level of this effective two- 
level system. The expression for the expectation value of 
the current in the PBCQ as well as the response of the 
tank circuit, weakly coupled to a PBCQ, have also been 
obtained. In Sec. Ill we describe the samples fabrication 
and the measurement setup. Comparison between theory 
and experimental data is discussed in the Sec. IV. 



II. THEORY 

A. Interaction of a phase-biased charge qubit with 
microwave power 

The phase-biased charge qubit (PBCQ), schematically 
shown in Fig. ^ consists of two Josephson junctions 
closed by a superconducting ring. The charge en of the 
island between the junctions is controlled by the gate 
voltage Vg via the capacitance C^, namely by the pa- 
rameter Ug = CgVg/e; eug is the polarization charge 
on the island. The junctions are characterized by the 
Josephson energies Eji, Ej2 and the phase differences ^i, 
82- The relevant energy values are the island's Coulomb 
energy, Ec = e^/2Ctot, where Ctot is the total capac- 
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itance of the island, and the effective Josephson en- 

1/2 

ergy sj = [Ej^ + ^ ^EjiEj2 cosS) . An impor- 
tant feature of the qubit is that its Josephson energy can 
be controhed by the external magnetic flux piercing 
the ring. In this paper, the ring inductance L is assumed 
to be small. Consequently, the total phase difference, 
S = Si -\- is approximately equal to = 27r<l>e/^o- 




FIG. 1: Scheme of the PBCQ 

The PBCQ is characterized by the Hamiltonian^: 

H = 4Ec{n - ng/2f - Eji cos Si - Ej2 cos^s, (1) 

which after quantization results in the following Hamil- 
tonian, written in the representation based on the eigen- 
states of the operator n, that is in the basis of the charge 
states |n): 



H = ^4£;c(n-n^/2)2|n)(n| + 



(2) 



A 



B 



^(|n + l)(n| + |n-l)(n|) + 



-^(|n + l)(n|-|n-l)(n|; 



A: 



s s 

-(Eji+Ej2)cos-, B = -{Eji-Ej2)sm-. (3) 



This Hamiltonian in the two-level approximation can be 
rewritten in the basis of the charge states {|1) , |0)]^: 
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(4) 



where the irrelevant term containing the unity matrix 
was omitted; 



C = 4Ec{l-ng), (5) 
andr^ are the Pauli matrices: |1) = |0) = — |0). 



We consider two possibilities for the excitation of a 
PBCQ: (a) via gate voltage: 

ng{t) = Ug -\-ng siucjt, S = const ^ (6) 

and (b) via magnetic flux: 

S{t) = ^0 + Ss'mut^ rig = const. (7) 

We shall consider flrst the time-independent case 
(which we denote by the subscript "0"). The eigenstates 
of the time-independent Hamiltonian 



Ho = -r. 



Co, 



-Tz, 



(8) 



denoted by {|— ) , |+)}, are related to the charge states 
{|0) , |1)} by the relation: 



I-) 
1+) 



S 



|0) 

ID 



Here 



S 



cos ^ 
-e-^^ sin I 



e^^ sin f 
cos ^ 



where the mixing angles r^, t/j are given by: 
sin T] = sj/ AE, cos r] = Co/ AE, 



sinV^ Bo/sj, cosV^ -Aq/ej 



with 



ej = ^/Al + Bl 



(9) 



(10) 



(11) 



(12) 



(13) 



(14) 



[4^c(l 



E% + 2EjiEj2 cos(5o. 



The diagonalization results in a Hamiltonian in the eigen- 
state basis: 



(15) 



where we denote the Pauli matrices, which operate in 
the eigenstate basis, by a^, so that we have |+) = 

In order to get the probabilities of the system to be 
in the eigenstates of the stationary Hamiltonian Hq, we 
rewrite the time- dependent Hamiltonian H{t) in this ba- 
sis. This Hamiltonian will be used in Sec. IV to solve 
the Bloch-type equation (the master equation) for the 
density matrix, whose diagonal elements deflne these 
probabilities^^. 
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We can split the Hamiltonian H{t) into two parts: 

%)=^o + ^i(t), (16) 

which gives: 

A P 

H\t) = S-^H{t)S = —a, + S-^Hi{t)S. (17) 

We consider first case (a), where the gate voltage is the 
time-dependent parameter as in Eq.Q. From Eq. (|T7)) 
it follows: 



^a(^) — - "^EcUg sin out X 



(18) 



X [cos T] -dz — sin r] sin ip -a^ -\- sin 77 cos ip • a^^] . 

We consider now case (b), where the magnetic flux is 
the time-dependent parameter, as in Eq. {TJ. In this 
case the time-dependent Hamiltonian ([TB|) can be rewrit- 
ten by expanding the quantities A and 5 in a Fourier 
series. However, we are interested only in small time- 
dependent perturbations and for the description of the 
experimental results we restrict ourselves to this case, 
when (5o=7r, ^/2<C7r. In the second approximation in 
5 the Hamiltonian can be written as 

Kit) = ^a, + (^ji+^j2)^sincjt-a,+ (19) 



B. Qubit-tank circuit arrangement 



/y. For our system we have&S: / = —Iq^z^ where the 
ground state current is 



e sin S 



(21) 



Thus, / = IqZ, where Z = (a^) = Sp(pdz)^ and p is 
the reduced density matrix. This means that the current 
flows with the probability P+ in one direction and with 
the probability P_ = 1 — P+ in the other direction; and 
the introduced value Z is equal to Z = 1 — 2P+. Thus 
from Eq. (|^nj) for the time-averaged phase shift a we 
have: 



2e 



tana k^Q • — L ■ { ^ - Z ^ Iq - ^ 



h 



dZ 



85 



85 



(22) 



where the bar stands for the time- averaging. We note 
that for a weakly driven system the function P+(cc^, ^, Ug) 
has the maxima (resonant peaks) at AE{5, Ug) = K ■ hco 
{K is an integer) and, consequently, its derivative (see 
Eq. (|22|)) has hyperbolic-like behavior. 

In the particular case, when ^ = tt we have Iq = and 
from Eq. it follows: 



tana : 



-A 



X = k' 



2e^LEc EjiEj2 



(23) 



£2 



(24) 



This means that the dependence of a on Ug at 5 
contains resonances at AE{ng) = K • hw. 



The current in a qubit ring is probed by the tank cir- 
cuit, which is weakly coupled through a mutual induc- 
tance M to the PBCQ^^-^"^. The PBCQ is characterized 
by the inductance LpscQ = L -\- Lj^ where L is the 
ring's inductance and Lj is an inductance defined by 
Lj^ = (2e/h) 81 / 85. The effect of this inductance on 
the tank circuit can be represented by an effective induc- 
tance: Lt Leff = Lt -\- M'^ / LpBCQ (where M is the 
mutual inductance) 

The experimentally measurable value is the phase shift 
between the voltage and current in the tank circuit a. 
The expression for the phase shift at the resonant fre- 
quency uot = 1/VLtCt is (see e.g. in Ref. [26]): 

t8inac^k^Q (20) 
Lj 

Here Q'^ = ujtCtRt, k^ = M^/{L • Lt), and we 
have neglected the ring's inductance L in the denom- 
inator. Thus the phase shift a is defined by Lj, i.e. 
by the current-phase relation I (5). This can be used 
to define the current-phase dependence in Josephson 
junctionsSSiSi. In the quantum case the current is equal 
to the expectation value of the current operator: / = 



III. SAMPLE PREPARATION AND 
MEASUREMENTS 

A SEM image of the gradiometer-type charge qubit 's 
core with two closely spaced mesoscopic Josephson junc- 
tions is shown along with the electrical circuitry in Fig.Q] 
The junction areas are slightly larger than 100x100 nm^ 
leading to critical currents Ic of 50-100 nA, which can 
be estimated from the measured tunnel resistances. The 
Cooper-box charge (on the island between the Josephson 
junctions) can be continuously varied by the gate volt- 
age. Both the junctions and the island are imbedded in a 
macroscopic (O.Sxlmm^) superconducting gradiometer- 
type- loop, which was done in order to minimize vibra- 
tion and magnetic noise. The single-Cooper-pair transis- 
tor and the loop were fabricated by e-beam lithography 
and shadow evaporation of aluminum. One loop of the 
gradiometer is inductively coupled by a flip-chip config- 
uration to the niobium high quality tank coil. 

We study the "qubit + tank" impedance as a function 
of the polarization charge eug and the phase difference 5^ 
by making use of a well-known impedance measurement 
technique. The tank circuit is driven by an rf current Irf 
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of frequency ljt close to the resonance frequency of the 
tank circuit. The phase difference a of the tank voltage 
(with respect to the phase of the applied current 7^/) is 
measured as a function of the gate voltage Vg = Vg{ng) 
and of the external magnetic flux = *^*e(^^)- These 
measurements show a shift in the resonant frequency of 
the "qubit + tank" arrangement due to a change in the 
effective inductance of the sample. The tank voltage was 
sequentially amplified by means of a cryogenic rf pream- 
plifier, a room temperature amplifier, and further rele- 
vant standard electronics. 

The measurements were carried out in a dilution-type 
refrigerator at a nominal temperature of about 10 mK. 
In order to minimize the noise level inherent in the to- 
tal gate voltage, we equipped the transistor's gate line 
(a ThermoCoax between 2 K and 10 mK) with conven- 
tional low-pass RC and microwave copper powder filters 
(Fig. Q). For an efficient thermalization of the charge 
gate, three microwave filters were mounted on different 
low-temperature plates (2 K, 50 mK, and 10 mK) of the 
refrigerator. The power attenuation of these 10 cm long 
filters was determined as a function of the frequency (up 
to 45 GHz) at room temperature. From the measure- 
ment results, we concluded that the total attenuation of 
this line was more than 80 dB in the GHz range. A high- 
attenuation ThermoCoax line along with two (on the 2 K 
and 10 mK) cooled commercial 20 dB-attenuators were 
used for applying a microwave power ( "UHF gate" ) to the 
sample. This microwave line is coupled inductively to the 
Pb-shield resonance cavity with the qubit inside. At the 
microwave frequency the current fed into the Pb shield is 
amplified by the quality factor of the resonator, produc- 
ing an electromagnetic field. In our measurements special 
care was taken to avoid the magnetic coupling between 
the microwave line and the qubit: (i) the qubit sample 
was placed across a Pb resonator for maximum electric 
(E)-field interaction and (ii) the gradiometer-type topol- 
ogy of the qubit circuit prevents the sample from the 
interaction via mutual inductance with the microwave 
line. As a further evidence for E-coupling via the length 
of the Al thin film charge-gate electrode, only the noise- 
like output signal on microwave power was obtained with 
sample after mechanical break of this thin film gate line. 

By passing a dc bias current through the tank coil 
(Fig. ^ we could simply control the flux-induced cur- 
rents circulating in the qubit ring, because of the mu- 
tual inductance between qubit and tank. For the tank, 
we prepared square-shaped Nb pancake coils on oxidized 
Si substrates^^. For flexibility, only the coil was made 
lithographically. We use an external capacitance Ct to 
be able to change the resonance frequency of the tank 
which, in this particular case was 28.9 MHz. The tank 
circuit was coupled by a 30 cm long piece of two- wire line 
to the cold HEMT- amplifier. Changes of the phase of the 
tank voltage oscillations due to variations of the Joseph- 
son inductance of the sample were measured by means of 
an averaging procedure: every measurement point was 
taken (with a time constant of 0.1 ms) 50 times and av- 



eraged. Cryogenic /i-metal and superconducting shields 
protected the sample against external magnetic and elec- 
tric noises. However, we could not take any special action 
to avoid the drift of electrostatic carriers within the sub- 
strate, so the 1/f noise due to background charge motion 
is not completely negligible in our experiments. 



IV. RESULTS AND DISCUSSION 

A. Analysis of the experimental results 

The dependence of the phase shift a on the gate volt- 
age is shown in Figs. I2K and In the first set of ex- 
periments, we used different frequencies of the microwave 
excitation at a nominally fixed power (see Fig.I2t)- In a 
similar manner to the results reported before^^ the Q^(n^) 
dependence exhibits clear peaks. Their positions depend 
on the frequency of excitation. Recently it was shown, 
that the peaks are due to resonant excitations of the sys- 
tem from ground to upper states^^. In the obtained de- 
pendencies a second set of the peaks is clearly seen (the 
grey arrows in Fig. These "additional" peaks would 
be due to two-photons excitation. In order to clarify this 
issue we fixed the frequency of the excitation and mea- 
sured a{ng) for different microwave powers. Indeed, as 
was expected, the "additional" peak structure becomes 
clearer for higher powers (see Fig.jSt^). Moreover, an ad- 
ditional structure appears in the a{^e) dependence as 
well (see Fig. |5Ji). 
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FIG. 2: Resonant excitation of the PBCQ: dependence of the 
phase shift a on the time-independent part of the dimension- 
less gate voltage at ^ = tt. The curves correspond to a 
fixed power of excitation in experiment (a) and an amplitude 
rig c^i 0.3 in theory (b). The varied parameter is the frequency 
a;/27r, which from the bottom to top curves is: 6.5, 7.1, 8.1, 
9.1 GHz. Upper curves are shifted vertically for clarity. Black 
(grey) arrows show the one- (two-) photon resonances. 

Let us extract the value of the minimum energy level 
separation A^min = Eji — Ej2 from the experimental 
results. In order to do that we define the position of 
the resonances, marked with the arrows in Fig. they 
correspond to AE{ng) = K - hu. We put these points 
in the Ug — A£^-plane: see Fig. ^ where the circles and 
squares correspond to one- and two- photon resonances, 
respectively, for which K = 1,2. The fitting of these data 
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FIG. 3: The same as in Fig. 2 but at a fixed frequency a;/27r 
8.2 GHz witfi tfie varied parameter: in experiment (a) being 
power of excitation (from bottom to top: -75, -63, -49, -42 
(dB)) and in theory (b) being ampHtude rig (from bottom to 
top: 0.1, 0.2, 0.3, 0.4). Upper curves are shifted. 



with the expression 



AE{ng,S = tt) 



[4^c (1 



{Ej^-Ej2Y 



(25) 

allows us to estimate both Eq and AiSmin = AE{ng = 
1, 5 = tt) = Eji - Ej2 (see Fig. ij. 




FIG. 4: Energy level separation as a function of Ug at 
6 = TT. Squares and circles correspond to one- and two- pho- 
ton resonances; solid line is a fit with Ec/h — 5 GHz and 
AEmin/h = {Eji - Ej2)/h = 5.5 GHz. 



B. Numerical calculations 

In this subsection we present the results of the quan- 
titative description of the system. In order to obtain the 
dependence of the tank voltage phase shift a on the sys- 
tem's parameters we made use of Eqs. (|22|) and (123- 
The probability that the upper level is occupied, P+(t), 
was obtained from the solution of the master equation 
for the density matrix as described in Ref.^^. In or- 
der to take into account the relaxation and dephasing 



processes the corresponding rates Treiax and F^^ 



are m- 
29. We 



eluded in the master equation phenomenologically 
note that the numerical solution of the master equation 
is the general approach used to describe the nonlinear 



dynamic behavior of a two-level system subjected to an 
external field of arbitrary amplitude and frequency (see 
in^^ and references therein). When the amplitude of the 
field is small and its frequency is close to the energy level 
separation devided by an integer, the analytical consid- 
eration, known as the rotating-wave approximation, can 
be applied to the description of the multiphoton tran- 
sitions (see e.g.^^). The latter approach was developed 
for the description of the PBCQ both analytically^^ and 
numerically^^ . 

In order to unify expressions (|T8j) and (|T9j) and to get 
the equations needed for numerical calculations, we write 
down the Hamiltonian H' as follows: 



rw 

H = -cr^ + -cr^ + -cr^. 



(26) 



Consequently, the evolution of the reduced density ma- 
trix p taken in the form 



^=2 



1 + Z X-iY 
X + iY l-Z 



(27) 



is described by the master equation in the form of Bloch 
equations (see inS^^): 



— = -Z 

dt h 

dY 

'dt 
dZ_ 

H 



-Y - r^x 



-X - r^r. 



s 



y-jrX- Trela. {Z " Z (0)) . 



(28) 
(29) 
(30) 



From these equations we get Z{t) which defines the oc- 
cupation probability of the upper level |+), P-\-{t) = 
P22{t) = ^(1 ~ ^(^))- We choose the initial condition 
to be X(0) = Y{0) = 0, Z(0) = 1, which corresponds to 
the system being in the ground state |— ). 

Quantitative analysis with Eqs. (|19|) and (|2HI)-(|Sn|) has 
shown that the case when the magnetic flux is time- 
dependent (Eq. O) is not consistent with the experi- 
mental results presented in this work. By that we have 
confirmed the argument, presented in Sec. Ill, that the 
qubit mainly is not excited via the magnetic flux, but 
rather via the gate voltage. So in what follows we will 
consider the case when the time-dependent parameter is 
the gate voltage (Eq. ©). Then the system is described 
by the Hamiltonian of Eq. (|18|) . 

We consider first the dependence of the phase shift a 
on Ug dit S = 7T. The results of the numerical calcu- 
lation are shown in Figs, (^b and In order to fit 
the experimental curves, by making use of Eqs. 
and we have taken: r^/{Ec/h) - 0.3 and 

^reiax/ {Ec /h) ~ 0.05 (which corresponds to the follow- 
ing decoherence and relaxation times: = F^^ :^ 0.7ns 

and Treiax = ^ relax — 4^^) and A = 0.1. This value of A 
is in good agreement with the value estimated from Eq. 

foi" the system's parameters experimentally accessi- 
ble. The relaxation and decoherence rates were assumed 
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to be independent of the system's parameters for simplic- 
ity. We note that the shape of the curves, in particular 
the widths and the heights of the resonances, is defined 
by three parameters: amplitude {rig) and the relaxation 
and decoherence rates. These values can be determined 
from the analysis of the widths and the heights of the 
resonances as e.g. in Ref.^^. But we rather fit the whole 
curves, which allows us to determine the system's param- 
eters. 

Now we consider the dependence of the phase shift a 
on 6 by making use of Eqs. and ((^ - ((5n|) : see Fig. 
El From the above considerations we have the following 
parameters: A, Ec and Eji — Ej2. But at (5 7^ tt we 
also need ^ji,2 (see Eq. (|^ ). At (5 7^ tt in the a — S- 
curve, due to the domination of the second term in Eq. 
(|22|) . the multiphoton resonances result in the hyperbolic- 
like behaviour with a at AE = Khw. We note 
that in the vicinity of (5 = tt the first term in Eq. 
decreases the value of a, which explains why the one- 
photon hyperbolic-like excitation is not symmetric about 
a = axis. From the position of these points, marked 
with the arrows in Fig. [5Ji (namely from the relation 
AE{ng,S) = Khuo, see Eq. (O), we found: Eji/h 
8Ec/h = 40 GHz, Ej2/h 6.9Ec/h = 34.5 GHz, and 
also Ug 0.85. With these values we have calculated 
the dependence of a on shown in Fig. [SJd. To fit 
the experimental curves we have taken: V^f)/ {Ec /h) = 
0.05 and T relax / {Ec /h) = 0.03, which correspond to the 
following decoherence and relaxation times: = F^^ 

AnS and Trelax = ^ relax " 




FIG. 5: Dependence of the tank voltage phase shift a on 
the phase difference 5. The curves correspond to the fixed 
frequency a;/27r = 7.05 GHz with the varied parameter: in 
experiment (a) being power of excitation (from bottom to 
top: -80, -60, -57 (dB)) and in theory (b) being amplitude 
rig (from bottom to top: 0.1, 0.2, 0.4). Upper curves are 
shifted. The arrows show the appearance of 1-, 2-, and 3- 
photon excitations at A^(^) = Khuo^ K = 1, 2, 3. 



V. CONCLUSIONS 

Multiphoton (namely, one-, two-, and three-photon) 
excitations of the PBCQ were observed experimentally 
and described theoretically. The multiphoton transitions 
manifest themselves in the dependence of the tank volt- 
age phase shift a on the qubit's parameters as follows: 
there are resonances in the dependence of a on Ug at 
6 = TT and there are hyperbolic-like dips and peaks in the 
dependence of a on 6. Theoretical fitting has allowed 
us to find out the qubit's parameters, particularly the 
relaxation and decoherence rates, which characterize the 
decoherence processes in the system. 
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